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Abstract 

CH4  adsorption  was  studied  experimentally  and  theoretieally  on  ZSM-5,  MOR,  and  ZSM- 
12  zeolites  using  calorimetric  measurements  at  195  K  and  plane  wave  DFT  calculations. 
Differential  heats  measured  on  four  different  H-ZSM-5  samples  were  determined  to  be  22.5±1 
kJ/mol,  independent  of  Bronsted  site  density  or  defect  concentration.  However,  DFT  calculations 
performed  using  various  functionals  and  on  the  most  stable  Bronsted  site  indicated  that  CH4  should 
bind  to  this  site  by  an  additional  1  to  7  kJ/mol,  a  discrepancy  that  is  due  to  the  inability  of  standard 
DFT  methods  to  capture  hydrogen-bonding  effects  accurately  with  CH4.  Differential  heats  for  CH4 
in  MOR  were  30±1  kJ/mol  at  low  coverages,  falling  to  25  kJ/mol  for  coverages  above  one 
molecule  per  8-membered-ring  side  pocket,  while  differential  heats  on  ZSM-12  were  initially  22.5 
kJ/mol,  decreasing  to  21  kJ/mol  with  coverage.  DFT  calculations  on  the  siliceous  form  of  the 
zeolites  were  able  to  predict  these  values  within  5  kJ/mol  in  most  cases.  The  results  indicate  that 
CH4  is  an  excellent  probe  molecule  for  characterizing  the  pore  structure  of  zeolites. 

Introduction 

Methane,  the  main  component  of  natural  gas,  would  be  an  attractive  fuel  for  vehicles  if  it 
could  be  stored  at  high  densities,  without  requiring  high  pressures  or  cryogenic  temperatures.  One 
method  for  achieving  high  storage  densities  at  reasonable  temperatures  and  pressures  involves 
adsorbing  methane  on  a  porous  solid  and  a  great  deal  of  effort  has  gone  into  developing  both 
materials  and  theoretical  understanding  for  this  application'’^.  Because  adsorption  of  methane  is 
physical  in  nature,  adsorption  isotherms  and  energetics  depend  in  a  rather  simple  manner  on  the 
composition  of  the  solid  adsorbent^  and  on  the  size  and  shape  of  the  pores  that  make  up  the  solid. 
Furthermore,  since  the  theory  for  adsorption  of  methane  is  reasonably  well  developed,  it  is  possible 
to  predict  what  the  adsorption  properties  of  a  crystalline  material  will  be  and  what  the  best  porous 
solid  structures  will  be"'. 

For  similar  reasons,  it  should  be  possible  to  use  experimental  data  to  characterize  the 
nanopores  in  unknown  or  poorly  crystalline  materials.  For  example,  the  zero-coverage  heats  of 
adsorption  for  methane  in  siliceous  zeolites  with  one-dimensional  pores  made  up  of  10-,  12-,  and 
14-membered  rings  (MR)  (TON,  MTW  (also  known  as  ZSM-12),  and  UTD-1)  were  reported  to 
be  27.2,  20.9,  and  14.2  kJ/mol,  respectively^.  If  one  knew  that  a  particular  siliceous  material  had 
cylindrical  pores,  the  heats  of  adsorption  could  be  used  to  estimate  the  pore  size.  Real  zeolite 
structures  have  more  complex  pore  structures  and  the  concept  of  using  the  distribution  of 
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adsorption  energies  eould  allow  one  to  map  the  strueture.  A  partieularly  interesting  material  for 
understanding  the  effeet  of  strueture  on  adsorption  is  the  zeolite  with  the  MOR  strueture.  MOR 
has  one-dimensional,  12-MR  ehannels  with  8  MR  side  poekets.  MOR  exhibits  unique  eatalytie 
properties  for  the  earbonylation  of  dimethyl  ether^,  most  likely  beeause  the  8-ring  poekets  are  able 
to  stabilize  the  intermediate  struetures. 

One  eomplieation  for  determining  pore  strueture  from  a  distribution  of  adsorption 
energies  is  that  most  zeolites  are  not  defeet  free  or  purely  silieeous.  Many  zeolite  struetures, 
ineluding  MOR,  are  synthesized  with  framework  A1  and  must  have  non-framework  eations  to 
balanee  the  eharge.  The  simplest  eharge-balaneing  eation  is  the  proton,  whieh  leads  to  Bronsted 
aeidity.  Even  simple,  polar  moleeules  (e.g.  CO)  will  adsorb  preferentially  on  these  sites  and  give 
an  added  eontribution  to  the  heats^.  With  non-polar  moleeules,  the  situation  is  less  elear. 
Interaetions  between  moleeular  O2  and  Bronsted  sites  are  reported  to  be  negligible^,  but  there  is 
both  experimental*’^  and  theoretieak*’’'^  evidenee  that  adsorption  on  the  Bronsted  sites 
eontributes  to  the  heats  of  adsorption  for  small  alkanes.  Beeause  adsorption  energies  have  been 
shown  to  seale  with  proton  affinities  for  both  protonated  and  hydrogen-bonded  adsorption 
eomplexes'^’^*  and  beeause  methane  has  a  mueh  lower  proton  affinity  eompared  to  even  ethane 
(543.5  kJ/mol  versus  596.3  kJ/moh"^,  the  influenee  of  aeid  sites  on  methane  adsorption  is  not  well 
known.  There  are  a  few  reported  studies  of  methane  adsorption  on  H-ZSM-5^^“^’,  but  the  role  of 
Bronsted  sites  in  adsorption  of  methane  has  not  been  elearly  demonstrated. 

In  order  to  measure  the  speetrum  of  adsorption  energies  in  a  material  that  has  more  than 
one  type  of  site,  it  is  essential  that  ealorimetrie  measurements  be  performed  at  an  appropriate 
temperature'*.  If  the  measurement  temperature  is  too  low  so  that  adsorption  is  not  reversible,  the 
adsorbate  will  not  be  able  to  sample  all  possible  sites.  In  this  ease,  if  adsorption  is  irreversible, 
eaeh  gas  dose  will  saturate  different  parts  of  the  sample,  leading  to  differential  heats  of  adsorption 
that  are  eonstant  with  eoverage,  even  if  the  sites  are  not  identieal'^.  On  the  other  hand,  for  weakly 
bound  moleeules  like  methane,  adsorption  should  be  performed  at  temperatures  that  are  low 
enough  for  adsorption  to  oeeur  preferentially  at  the  strongest  sites.  For  example,  if  a  material  has 
two  types  of  sites  (e.g.  in  MOR,  these  might  be  poeket  sites  and  main-ehannel  sites),  separated  in 
energy  by  values  similar  to  fcT  (where  k  is  the  Boltzmann  eonstant  and  fcT  -'2.5  kJ/mol  at  room 
temperature),  a  signifieant  fraetion  of  moleeules  added  to  a  fresh  sample  will  oeeupy  the  weaker 
sites  at  equilibrium,  espeeially  if  the  weaker  sites  are  present  in  higher  eoneentrations  and  if 
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adsorption  on  those  sites  is  entropieally  preferred.  Indeed,  sample  ealculations  have  suggested  that 
it  is  likely  impossible  to  observe  minority  sites  having  a  heat  of  adsorption  that  is  5  kJ/mol  higher 
than  that  of  the  majority  sites  using  differential  ealorimetry  at  room  temperature'^. 

Past  work  from  one  of  our  laboratories  demonstrated  that  site  heterogeneity  in 
physisorption  can  be  measured  using  differential  calorimetry  at  195  K^.  In  siliceous,  defect- free 
ZSM-5  (a  zeolite  with  the  MFI  structure),  the  differential  heats  of  adsorption  for  CO  were  found 
to  be  16±1  kJ/mol,  independent  of  coverage.  On  the  analogous  acidic  H-ZSM-5,  the  differential 
heats  for  CO  adsorption  were  10  kJ/mol  higher  at  low  coverages,  falling  to  the  same  value  as  on 
the  siliceous  ZSM-5  at  coverages  above  one  CO  per  Bronsted  site.  Molecules  such  as  O2,  which 
do  not  interact  with  Bronsted  sites,  showed  the  same  heats  of  adsorption  on  both  samples^®. 
Therefore,  the  coverage-dependent  differential  heats  provide  a  measure  of  the  site  energetics  for 
these  materials. 

In  the  present  work,  we  set  out  to  determine  whether  differential  calorimetry  can  be  used 
to  distinguish  structural  differences  in  similar  zeolites,  using  MOR  and  MTW  as  the  examples. 
Using  both  theoretical  calculations  and  experimental  measurements,  we  will  show  that  methane 
adsorption  at  195  K  appears  to  be  an  excellent  probe  for  determining  structural  differences.  By 
studying  CH4  adsorption  in  a  series  of  H-ZSM-5  samples,  we  show  that  the  presence  of  Bronsted 
sites  does  not  affect  the  heats  of  adsorption  for  CH4. 

Methods 

Experimental 

The  Tian-Calvet  calorimeter  used  in  this  study  was  home-built  and  has  been  described  in 
detail  elsewhere'^.  The  instrument  was  constructed  from  five,  2.54-cm  square,  thermal-flux  meters 
(International  Thermal  Instrument  Company,  Del  Mar,  CA)  placed  between  a  cubic,  pyrex  sample 
cell  and  a  large  A1  block.  The  1-g  zeolite  samples,  pressed  into  wafers  and  placed  at  the  bottom  of 
the  pyrex  cube,  were  covered  with  quartz  chips  in  order  to  prevent  heats  losses  due  to  radiation 
out  the  top.  The  sample  was  evacuated  using  a  mechanical  pump  and  exposed  to  gases  from  a 
calibrated  GC  sample  loop.  All  experiments  were  performed  at  195  K,  maintained  by  placing  the 
Al-block  heat  sink  in  a  styrofoam  container  with  dry  ice.  The  thermopiles  had  previously  been 
calibrated  by  passing  current  through  a  Pt  wire  placed  between  the  sample  cell  and  the  thermopiles. 
The  uncertainty  in  our  measurement  of  each  point  was  2%,  corresponding  to  '-0.5  kJ/mol; 
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however,  the  uneertainty  in  our  calibration  factor  is  larger,  so  that  the  uncertainty  in  the  absolute 
values  of  the  differential  heats  is  ~1  kJ/mol. 

The  zeolites  in  this  study,  along  with  some  of  their  properties,  are  listed  in  Table  1.  Two 
of  the  H-ZSM-5  and  the  H-MOR  samples  are  commercially  available  from  Zeolyst.  The  two  ZSM- 
5  samples,  prepared  at  the  University  of  Delaware,  used  fluoride-containing  media  to  achieve  low 
silanol  defect  concentrations^^  A  typical  synthesis  was  as  follows:  distilled  water  was  mixed  with 
Ludox  AS-40  silica  sol,  tetrapropylammonium  bromide  (TPABr),  NH4F,  and  A1(N03)3  •  9H2O. 
(A  purely  siliceous  ZSM-5  was  prepared  using  the  same  procedure,  but  without  the  addition  of 
A1(N03)3  •  9H2O.)  The  mixture  was  heated  at  450  K  for  seven  days  in  a  Teflon-lined  autoclave. 
After  filtering,  the  solid  was  calcined  at  823  K  in  flowing  dry  air  for  6  h.  These  samples  were 
noticeably  more  hydrophobic  and  2-propanamine  molecules  in  excess  of  twice  the  Bronsted-site 
densities  could  be  easily  evacuated  at  room  temperature,  whereas  the  excess  amine  molecules  in 
the  commercial  samples  could  only  be  removed  by  heating.  The  H-ZSM-12  sample  was  used  in  a 
previous  study  and  has  been  described  there  in  more  detaifi.  Bronsted-site  concentrations  were 
determined  using  simultaneous  temperature-programmed  desorption  and  thermogravimetric 
analysis  of  2-propanamine^^.  Pore  volumes  were  determined  from  gravimetric  uptakes  of  n-hexane 
at  room  temperature  and  10  torr,  assuming  the  pore  volume  becomes  filled  with  liquid-like  n- 
hexane. 

Calculations 

All  electronic  structure  calculations  were  carried  out  using  Vienna  ah  initio  simulation 
package  (VASP)^^’^"^-  a  plane  wave  density  functional  theory  (DFT)  code.  PAW^^  potentials  with 
generalized  gradient  approximation  (GGA)  using  the  PBE  exchange  correlation  functionafi^  with 
Grimme-D2^^  semi-empirical  dispersion  corrections  were  used.  The  Brillouin  zone  was  sampled 
using  only  the  gamma  point  in  view  of  the  large  supercell  size.  A  Gaussian  smearing  of  0.1  eV 
was  used  and  all  energies  were  extrapolated  to  0  K.  A  plane  wave  cutoff  of  400  eV  and  density 
wave  cutoff  of  600  eV  were  used  with  a  convergence  criteria  of  all  forces  being  smaller  than  0.02 
eV/A.  We  refer  to  this  level  of  theory  as  DFT-D2  in  this  work.  The  calculated  (experimental^  in 
parenthesis)  lattice  parameters  for  the  three  zeolites  in  their  siliceous  forms  are:  (a)  20.2  x  19.9  x 
13.3  A  (20.1  X  19.7  X  13.1  A)  for  the  orthorhombic  phase  of  ZSM-5,  (b)  18.4  x  20.7  x  7.6  A  (18.3 
X  20.5  X  7.5  A)  for  the  orthorhombic  mordenite  (MOR)  structure,  and  (c)  25.8  x  5.3  x  12.2  A  (25.6 
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X  5.3  X  12.1  A  with  P=109.3°  for  both)  for  zeolite  ZSM-12  in  the  monoelinic  form.  In  these 
ealeulations,  the  atoms  of  the  zeolite  were  fixed  at  eoordinates  speeified  in  the  structural 
information  published  by  IZA  website  All  molecular  adsorption  calculations  for  ZSM-5  and  MOR 
were  performed  in  a  supercell  comprising  one  full  periodic  unit  while  those  on  ZSM-12  were 
carried  out  in  a  supercell  that  had  two  unit  cells  (two  units  in  the  y-axis);  these  cells  enable  us  to 
explore  low  adsorbate  coverages.  All  atoms  were  relaxed  in  these  calculations  and  spin 
polarization  was  considered  for  the  oxygen  molecule  in  the  gas  and  adsorbed  phase.  Gas  phase 
calculations  were  carried  out  in  the  same  level  of  theory  as  above  in  a  cell  of  size  20  x  20.5  x  21 
A  allowing  for  more  than  10  A  of  vacuum  between  periodic  images.  Figures  of  zeolite  structures 
and  adsorbate  configurations  were  generated  using  VESTA^^. 

The  binding  energy  (BE)  of  an  adsorbate  molecule  is  calculated  as: 

~  ^ molecule+Zeolite  ^molecule,  gas  ^Zeolite  (1) 

where  ^’^ofecuie+zeofite  energy  of  the  zeolite  with  the  adsorbed  molecule, 

^molecule  gas  energy  of  the  gaseous  molecule,  and  Ezeoiue  is  the  total  energy  of  the 

zeolite  (catalyst).  We  report  results  on  both  Bronsted  and  non-Bronsted  forms  of  multiple  sites  on 
the  zeolite  ZSM-5;  for  MOR  and  ZSM-12,  only  the  non-Bronsted  sites  are  considered.  The 
calculation  of  and  on  a  Bronsted  site  involved  performing  electronic 

structure  calculations  on  a  zeolite  supercell  with  a  single  Si  atom  being  replaced  by  an  A1  atom 
(and  adding  one  H  atom  to  a  neighboring  oxygen  atom);  for  the  corresponding  non-Bronsted  site, 
similar  calculations  were  performed  on  its  topologically  identical  site  in  a  siliceous  supercell. 

Zero  point  energy  (ZPE)  corrections  were  included  (unless  otherwise  stated)  in  all  binding 
energy  values  reported  here.  The  vibrational  frequencies  Vj  of  the  adsorbate  (and  the  Bronsted 
proton,  if  present)  were  computed  using  harmonic  approximation  and  ZPE  was  calculated  using 
the  following  formula 

=  (2) 
Since  there  are  multiple  topologically  identical  locations  in  a  zeolite  supercell,  we  also 
performed  sample  calculations  to  test  if  the  presence  of  a  Bronsted  site  in  a  supercell  affects 
adsorption  on  a  topologically  identical  non-Bronsted  site  and  compared  this  with  the  binding 
energy  values  on  the  siliceous  supercell;  we  observed  negligible  (<  1  kJ/mol)  differences  in  all 


cases. 
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Temperature  eorreetions  were  added  to  the  binding  energy  to  get  enthalpy  of  adsorption. 
This  was  done  by  calculating  temperature-dependent  entropy  values  for  the  adsorbate  in  the  gas 
and  adsorbed  phases,  fitting  a  polynomial  (Shomate  equations)  to  this  data,  calculating  heat 
capacity  from  the  resulting  expression,  and  calculating  a  temperature  contribution  to  the  total 
energy.  The  details  of  this  procedure  are  given  in  the  supporting  information  (Method  SI). 

Results  and  Discussion 
H-ZSM-5 

Experimental 

In  order  to  understand  the  possible  effects  of  Bronsted  sites  and  silanol  defects  on  CH4 
adsorption,  calorimetric  measurements  were  performed  at  195  K  on  the  four  H-ZSM-5  samples 
listed  in  Table  1.  The  differential  heats  are  reported  as  a  function  of  coverage  in  Figure  1,  with 
results  for  the  two  Zeolyst  samples  in  Figure  la)  and  for  the  two  defect- free  samples  in  Figure  lb). 
Within  experimental  uncertainty,  the  heats  of  adsorption  at  coverages  below  0.6  mmol/g  were 
independent  of  coverage  on  all  of  the  samples.  The  differential  heats  were  22.5±1  kJ/mol  on  the 
two  Zeolyst  samples  and  22±1  on  the  two  defect- free  samples.  The  Bronsted-site  densities  on  these 
samples  ranged  from  0.010  mmol/g  to  0.47  mmol/g,  implying  that  interactions  with  the  Bronsted 
sites  were  too  small  to  observe  in  calorimetry.  The  very  small  difference  between  the  results  on 
the  defect-free  samples  and  the  Zeolyst  samples  were  less  than  the  experimental  uncertainty. 

Interestingly,  the  equilibrium  constants,  determined  from  the  slope  of  the  isotherm  near 
zero  coverage  and  reported  in  Table  2,  were  not  identical  on  the  four  samples.  The  values  for  the 
four  samples  differed  by  a  factor  of  about  4,  with  H-ZSM-5(50)  showing  the  highest  value  and 
ZSM-5(F,Si)  the  lowest.  There  is  a  strong  correlation  between  the  equilibrium  constant  and  the 
Bronsted-site  density  but  there  does  not  appear  to  be  any  dependence  on  whether  or  not  the  samples 
are  defect-free.  If  the  change  in  equilibrium  constant  were  simply  due  to  enthalpic  interactions 
with  the  Bronsted  sites,  the  factor  of  4  at  195  K  would  imply  that  the  heats  of  adsorption  should 
be  2.2  kJ/mol  higher  on  H-ZSM-5(50)  compared  to  ZSM-5(F,Si).  This  is  significantly  higher  than 
the  uncertainty  in  the  calorimetric  measurements,  implying  that  there  may  be  some  entropic  effects 
associated  with  adsorption  near  the  Bronsted  sites. 

For  example,  since  adsorbed  CH4  is  likely  very  mobile,  we  suggest  that  the  effect  of  the 
Bronsted  sites  might  be  to  reduce  this  mobility.  The  change  in  adsorption  entropy  associated  with 
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a  factor  of  4  change  in  equilibrium  eonstant  is  11.5  J/mofK  at  195  K.  If  Bronsted  sites  were  to 
reduee  the  mobility  from  something  like  a  two-dimensional  to  a  one-dimensional  gas,  the  partition 
funetions  eould  be  ealeulated  using  Equations  (3)  and  (4): 

5ira„.2D=fi[ln(5^)+in(“)  +  2]  (3) 

S,„„.1D  =  «  (log  (^3^)  +  log(L„)  +  1.  s)  (4) 

SA 

where  m  is  the  mass  of  the  moleeule,  k  is  the  Boltzmann  eonstant,  and  h  Planek’s  eonstant.  —  is 

the  average  surfaee  area  available  for  a  moleeule  at  the  referenee  state  of  a  monolayer  eoverage, 
equal  to  200  pm  x  600  pm  for  ZSM-5'\  and  Lq  can  be  taken  to  be  a  lattiee  parameter,  on  the  order 
of  600  pm  in  this  ease.  For  CH4  at  195  K,  the  differenee  between  these  values  is  nearly  20  J/mohK. 


Calculations 

Sinee  previous  theoretieal  studies  have  indieated  hydroearbons  should  interaet  with 
Bronsted  sites we  performed  theoretieal  studies  on  H-ZSM-5  to  understand  what  effeet 
Bronsted  sites  might  have  on  the  adsorption  of  simple  moleeules.  In  addition  to  examining  the 
adsorption  of  CH4,  we  also  examined  adsorption  of  CO  on  H-ZSM-5,  sinee  previous  measurements 
indieated  that  CO  adsorbs  on  Bronsted  sites  with  an  additional  10  kJ/mol  (This  result  was  repeated 
here  for  the  H-ZSM-5(F,A1)  sample,  with  data  shown  in  Figure  Sl)^. 

Two  loeations  for  the  aeid  sites  were  eonsidered.  The  most  stable  aeid  site  (A17-017-Si8 
defined  per  the  seheme  adopted  by  IZA)  is  shown  in  Figure  S2,  where  the  Bronsted  proton  (bonded 
to  017)  is  hydrogen  bonded  to  another  oxygen  atom  (0-H  distanee  of  1.85  A).  A  less  stable  site 
(Si 12-026-All  2)  (by  0.23  eV  or  -'22  kJ/mol)  was  also  eonsidered  at  the  entranee  of  the  sinusoidal 
ehannel  and  pointing  into  the  straight  ehannel  is  shown  in  Figure  S3.  These  two  sites  have  been 
eonsidered  in  the  literature  as  the  likely  positions  for  the  Bronsted  sites  in  HZSM-5^^“^"^.  The 
enthalpy  of  adsorption  at  the  DFT-D2  level  of  theory  for  CO  and  CH4  on  these  two  loeations,  with 
and  without  the  Bronsted  site,  are  reported  in  Table  3  and  diseussed  below.  The  eorresponding 
binding  energy  values  (ZPE  eorreeted)  are  in  Table  SI.  Figures  S4-S7  show  the  adsorption 
struetures  for  CH4  and  CO  on  the  Bronsted  and  non-Bronsted  forms  of  the  two  sites  eonsidered. 

CO  binds  via  the  earbon  atom  on  both  Bronsted  and  non-Bronsted  sites  of  ZSM-5  similar 
to  the  observations  made  in  eomputational  studies  of  other  proton-form  zeolites^^’^^;  in  the 
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presence  of  a  Bronsted  site,  the  molecule  points  towards  the  Bronsted  proton,  leading  to  an 
additional  stabilization  of  about  7  kJ/mol  for  T7-017-T8  and  22  kJ/mol  eV  for  the  T12-026-T12 
site.  The  adsorption  enthalpy  values  (-26.7  kJ/mol)  on  the  most  stable  Bronsted  acid  site  (T7-017- 
T18)  match  remarkably  well  with  the  microcalorimetric  data  (26.5  kJ/mol)  for  carbon  monoxide 
on  ZSM-S^'*.  On  the  other  hand,  the  binding  energy  on  the  Bronsted  site  of  T12-026-T12  is  about 
14  kJ/mol  higher  than  that  of  T7-017-T8. 

The  calculated  adsorption  enthalpies  of  CH4  on  the  Bronsted  (-23.2  kJ/mol)  and  siliceous 
forms  -21 .2  kJ/mol)  of  the  T7-017-T8  site  are  close  to  experimental  values  (~22  kJ/mol);  however, 
DFT-D2  calculations  show  an  effect  of  the  Bronsted  site  of  about  2  kJ/mol.  The  binding  energy 
on  the  Bronsted  site  of  T12-026-T12  is  about  12  kJ/mol  larger  than  that  on  T7-017-T8,  leading 
again  to  a  more  pronounced  effect  of  the  Bronsted  site  at  the  T12-026-T12  location.  It  should  be 
noted  that  the  lower  effect  of  hydrogen  bonding  for  the  T7-017-T8  site  is  potentially  due  to 
stabilization  of  the  Bronsted  hydrogen  atom  through  hydrogen  bonding  (Figure  S2).  Recently, 
Tuma  and  Sauer  adopted  a  hybrid  MP2:  DFT  method  to  calculate  the  binding  energy  and  proton 
transfer  barriers  of  methane  on  a  number  of  proton-form  zeolites^^.  This  higher  level  of  theory  can 
be  used  to  benchmark  our  results.  The  authors  report  a  ZPE-corrected  binding  energy  value  of  - 
28.5  kJ/mol  on  the  T7-017  site  (although  named  T7-07  according  to  their  naming  convention) 
compared  to  our  -26  kJ/mol  (see  Table  SI).  In  addition,  we  performed  additional  calculations  on 
the  T12  site  studied  by  Tuma  and  Sauer  (T 12-020  in  their  naming  convention  and  T 12-08  in 
ours);  the  ZPE-correted  binding  energy  values  of  methane  adsorption  on  this  site  was  found  to  be 
-29.9  kJ/mol,  in  close  agreement  with  the  published  value  of  -30.1  kJ/mol  (see  Table  S2).  The 
effect  of  Bronsted  site  for  methane  adsorption  on  the  T 12-08  site  is  about  3.5  kJ/mol. 

Because  experimental  data  suggests  acid  sites  have  a  negligible  effect  on  CH4  adsorption 
enthalpies,  we  investigated  the  origin  of  this  DET-predicted  effect  by  considering  other  exchange 
correlation  functionals,  as  well  as  alternative  treatments  of  dispersion  effects.  Table  4  lists  the 
different  levels  of  theory  considered  in  this  study,  and  the  corresponding  calculated  difference  in 
the  binding  energy  values  (not  ZPE  corrected)  of  CH4  on  the  Bronsted  and  non-Bronsted  forms  of 
the  T7-017-T8  site.  The  difference  values  are  between  0.01  and  0.07  eV  ('-I  to  7  kJ/mol)  in  all 
cases  and  arise  primarily  due  to  the  choice  of  the  exchange  correlation  functional  used  in  the 
calculations.  Similar  observations  can  be  made  for  the  T12-026-T12  site  as  well  (see  Table  S3). 
These  results  are  consistent  with  other  theoretical  studies^^’^®  but  do  not  agree  with  the 
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experimental  observations.  We  suggest  the  diserepancy  between  theory  and  experiment  in  this  ease 
is  due  to  the  inability  of  standard  DFT  methods  to  eapture  hydrogen-bonding  effeets  aceurately. 
Various  benehmarking  studies  on  the  aeeuraey  of  DFT  funetionals  to  treat  hydrogen  bonding  have 
been  reported"^"^^’.  While  the  specifie  conclusions  of  these  studies  vary  (some  report  under¬ 
prediction  while  others  observed  over-prediction  with  respect  to  a  higher  level  of  theory  such  as 
coupled  cluster  calculations)  and  the  reference  set  of  molecules  considered  are  different,  they  all 
point  to  intrinsic  errors  in  standard  generalized  gradient-corrected  exchange  correlation  functionals 
in  treating  hydrogen  bonding,  in  line  with  our  observations.  In  view  of  our  experimental  data,  the 
Bronsted  sites  will  not  be  considered  further  in  our  calculations. 

MORandMTW 

Experimental  Data 

In  order  to  determine  the  effect  of  zeolite  structure  on  adsorption  properties,  we  performed 
calorimetric  measurements  of  CH4,  O2,  and  Ar  on  MOR  and  ZSM-12.  As  discussed  in  the 
Introduction,  the  main  channels  in  both  MOR  and  ZSM-12  are  12  MR.  However,  MOR  also  has 
8-MR  pockets  on  the  sides.  Although  the  volume  of  these  pockets  is  small  compared  to  that  of  the 
main  channels,  the  corrugation  that  results  from  the  openings  to  the  side  pockets  causes  the  pore 
surfaces  to  be  rougher  on  the  molecular  scale.  As  with  CH4,  O2  and  Ar  are  not  expected  to  interact 
with  Bronsted  sites,  so  that  most  of  the  interactions  should  be  due  to  dispersion  forces  with  the 
siliceous  walls^'*. 

Figure  2  shows  the  differential  heats  for  CH4  on  ZSM-12  and  MOR.  On  ZSM-12,  the 
differential  heats  are  reasonably  constant,  starting  at  23.0  kJ/mol  at  zero  coverage  and  decreasing 
slightly  to  21.5  kJ/mol  at  0.4  mmol/g.  In  a  previous  calorimetric  study  carried  out  at  room 
temperature,  the  differential  heats  also  decreased  slightly  with  coverage  but  the  reported  heats  were 
about  2  kJ/mol  lower^.  About  1.2  kJ/mol  of  this  difference  can  be  explained  by  the  heat  of 
adsorption  changing  with  temperature.  This  value  was  determined  by  assuming  CH4  loses  its 
translational  degrees  of  freedom  in  going  to  the  adsorbed  phase,  so  that  the  heat  capacities  of  gas- 
phase  and  adsorbed-phase  CH4  change  by  (3/2)R.  In  any  case,  the  numbers  here  are  probably  more 
accurate  than  those  reported  previously  due  to  the  fact  that  lower  gas  pressures  were  required  to 
obtain  a  given  coverage  at  195  K.  Although  ZSM-12  consists  of  12-MR  channels,  the  differential 
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heats  are  very  close  to  that  found  on  the  10-MR  zeolite,  ZSM-5.  This  is  due  to  the  fact  that  the 
intersecting  channels  in  ZSM-5  make  the  channels  effectively  larger. 

The  heats  of  adsorption  for  CH4  on  MOR  were  significantly  higher  than  those  on  MTW. 
At  low  coverages,  the  differential  heats  in  MOR  were  almost  30  kJ/mol.  For  coverages  between 
0.25  and  0.45  mmol/g,  the  differential  heats  fell  to  25  kJ/mol,  then  remained  constant  at  that  value. 
The  higher  heats  are  almost  certainly  associated  with  the  8-MR  side  pockets.  While  there  are  no 
simple,  one-dimensional  zeolites  with  8-MR  channels,  it  is  noteworthy  that  the  zero-coverage 
differential  heats  were  only  about  27  kJ/mol  on  the  10-MR,  one-dimensional  zeolite,  TON. 
Furthermore,  the  coverage  at  which  the  differential  heats  fall  on  MOR  is  close  to  the  side-pocket 
concentration  of  0.35  mmol/g.  The  fact  that  there  is  no  sharp  demarcation  at  0.35  mmoFg  is  due 
to  the  adsorption  energies  for  molecules  in  the  side  pockets  and  in  the  corrugated  main  channels 
differing  by  a  relatively  small  amount,  so  that  both  sites  are  being  occupied  simultaneously  at  195 
K,  as  discussed  in  the  Introduction'^. 

Figure  3  shows  the  corresponding  differential  heats  for  O2  in  the  ZSM-12  and  MOR 
samples.  Heats  of  adsorption  for  O2  in  ZSM-12  were  16  kJ/mol,  independent  of  coverage.  As  with 
CH4,  this  value  is  again  similar  to  the  heat  of  adsorption  reported  for  O2  in  siliceous  ZSM-5^°.  The 
differential  heats  for  O2  on  MOR,  20  kJ/mol,  were  higher  than  those  found  on  ZSM-12;  however, 
unlike  the  case  for  CH4,  the  differential  heats  for  O2  were  independent  of  coverage.  Because  O2 
should  be  able  to  fit  in  the  8-ring  side  pockets,  we  suggest  that  difference  in  adsorption  enthalpies 
between  the  side  pockets  and  the  main  channels  is  similar  or  less  than  fcT,  so  that  both  sites  fill 
simultaneously.  The  results  for  Ar  in  MOR  and  ZSM-12,  shown  in  Figure  4,  are  nearly  the  same 
as  those  for  O2. 


Calculations 

The  adsorption  of  CH4,  O2,  and  Ar  in  siliceous  MOR  and  ZSM-12  was  studied  using  DFT- 
D2.  Figure  S8  shows  a  schematic  of  MOR.  Molecules  can  adsorb  in  the  12-MR  channel  and  the 
8-MR  side  pocket.  While  MOR  also  has  an  8-MR  channel  as  marked  in  Figure  S8,  the  dimension 
of  its  minor  axis  (4.3  A)  prevents  adsorption  in  this  region;  we  did  not  observe  strong  adsorption 
in  the  8  MR  channels  except  for  a  few  less  stable  structures  at  the  intersection  of  8  MR  channels 
and  8  MR  side  pockets.  The  12  MR  channel  close  to  the  walls  (marked  region  II)  and  the  region 
of  8  MR  side  pockets  (marked  region  II)  are  potential  locations  for  adsorption. 
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The  calculated  adsorption  enthalpy  of  CH4,  O2,  and  Ar  in  the  two  regions  of  MOR  are 
given  in  Table  3  (ZPE  corrected  binding  energies  are  in  Table  S4)  and  their  adsorption  structures 
are  shown  in  Figure  5.  CH4  binding  energies  vary  from  -28  kJ/mol  to  -17  kJ/mol)  while  O2  and 
Ar  binding  energy  variations  are  smaller.  In  all  cases,  DFT-D2  predicts  a  higher  binding  energy 
for  the  side  pockets  than  the  straight  channels.  Further,  the  DFT-D2  binding  energy  in  the  8-MR 
side  pockets  match  experimental  heats  at  low  coverage,  corroborating  our  argument  that  the  side 
pockets  get  filled  first.  DFT-D2  under-predicts  the  binding  energy  of  CH4  in  the  12  MR  channel 
by  5  kJ/mol  (compared  to  experimental  heats  at  high  coverage).  This  under-prediction  appears  to 
be  uniform  for  all  three  adsorbates  because  accounting  for  this  brings  DFT-D2  predictions  to 
within  2  kJ/mol  of  experimental  heats  at  high  coverages  for  O2  and  Ar.  This  uniform  deviation 
probably  arises  from  intrinsic  error  in  the  variation  of  the  dispersion  potential  as  a  function  of 
distance  of  the  adsorbate  from  the  wall.  The  trends  from  DFT-D2,  however,  suggest  that  the  heats 
of  adsorption  of  O2  and  Ar  are  independent  of  the  coverage  because  the  difference  between  the 
adsorption  enthalpies  of  the  side  pockets  and  the  main  channel  is  small  compared  to  fcT. 

A  schematic  of  ZSM-12  showing  the  12-MR  channels  is  given  in  Figure  S9.  The  DFT-D2 
adsorption  enthalpy  of  CH4,  O2,  and  Ar  in  ZSM-12  is  given  in  Table  4  (ZPF  corrected  binding 
energy  values  are  in  Table  S5)  and  the  adsorption  structures  are  shown  in  Figure  6.  CH4  (-20.3 
kJ/mol)  binds  more  strongly  than  O2  (-13.8  kJ/mol)  and  Ar  (-16.6  kJ/mol),  similar  to  what  was 
found  in  MOR;  these  are  within  4  kJ/mol  of  the  experimental  differential  heats.  The  binding  energy 
of  CH4  in  ZSM-12  is  within  1  kJ/mol  of  that  onZSM-5  in  agreement  with  the  experiments.  Clearly, 
comparing  the  three  zeolites,  the  binding  energy  of  adsorbates  is  higher  in  confined  spaces  of  the 
zeolite  (such  as  the  side  pockets  of  MOR)  due  to  increased  dispersion  interactions.  Overall,  the 
calculations  confirm  our  initial  proposition  and  experimental  results  that  a  probe  molecule  that 
interacts  with  the  zeolite  purely  through  dispersion  interactions  (such  as  methane)  can  be  used  to 
characterize  the  pore  geometry  of  zeolites. 


Conclusions 

Experimental  differential  heats  of  adsorption  for  CH4  depend  strongly  on  the  solid 
nanostructure  but  not  on  the  presence  of  Bronsted  sites,  implying  that  calorimetric  measurements 
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at  195  K  could  be  used  to  probe  the  nanostrueture  of  unknown  materials.  The  eomparison  of  the 
experimental  results  with  DFT  demonstrated  that  the  theory  aecurately  deseribes  the  effect  of 
zeolite  strueture  on  CH4  adsorption  but  over-prediets  the  effect  of  Bronsted  sites. 
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Table  1.  Zeolite  samples  used  in  this  study. 


Zeolite 

Si/Ab 

Bronsted-acid  site 
densities  (pmol/g) 

Pore  volume 
with  n-hexane 
(cmVg) 

Zeolyst,  H-ZSM-5(50) 

50 

470 

0.1856 

Zeolyst,  H-ZSM-5(280) 

280 

80 

0.1875 

Defect  free,  H-ZSM- 
5(F,A1) 

- 

240 

0.1744 

Defect  free,  ZSM-5(F,Si) 

- 

10 

0.1696 

H-MOR 

20 

- 

0.1026 

ZSM-12 

- 

20 

0.0939 

Table  2.  Equilibrium  constants  for  adsorption  of  CH4  on  the  four  H-ZSM-5  samples  at  195  K, 
determined  from  the  slope  of  the  isotherms  near  zero  coverage. 


Zeolite 

Slope  at  zero  coverage 
10‘^  (mol/(g*Pa)) 

H-ZSM-5(50) 

6.47 

H-ZSM-5(280) 

4.15 

H-ZSM-5(F,A1) 

4.80 

ZSM-5(F,Si) 

1.54 

17 
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Table  3:  Calculated  enthalpy  of  adsorption  for  CO  and  CH4  (in  kJ/mol)  on  Bronsted  and  non- 
Bronsted  forms  of  two  sites  in  ZSM-5  at  195  K. 


Adsorbate 

T7-017-T8‘‘ 

T12-026-T12^ 

Bronsted’’  Non- 

Bronsted’’  Non- 

Bronsted’’ 

Bronsted’’ 

CO 

-26.7  -19.1 

-41.7  -19.6 

Methane 

-23.2  -21.2 

-35.0  -20.1 

“  “T”  refers  to  the  tetrahedral  atom  that  is  either  Si  or  A1  atom;  ’’  “Bronsted”  refers  to  adsorption 
on  a  Bronsted  site  formed  by  replacing  a  Si  atom  of  Tx-Oy-Tz  site  by  A1  atoms  and  adding  a 
hydrogen  to  the  oxygen  atom  “y”.  For  T7-017-T8,  A1  replaces  a  Si  atom  in  the  tetrahedral 
position  7  and  for  T12-024-T12  A1  replaces  Si  atom  in  the  tetrahedral  position  12. 

Table  4:  Difference  in  CH4  binding  energies  between  Bronsted  and  non-Bronsted  T7-017-T8 
site  of  ZSM-5  using  different  functionals  and  dispersion  treatments 

Method 

Dispersion-corrected 

Difference  [eV 

(kJ/mol)]‘’’’’ 

Reference  (PBE+D2)^^'^^  Yes 

0.05  (4.7) 

PBE^’’ 

No 

0.05  (4.7) 

No 

0.05  (4.9) 

RPBE^o 

No 

0.01  (1.3) 

PBE^’’  +  DS"” 

Yes 

0.07  (6.6) 

PBEsol^^ 

No 

0.07  (6.5) 

BEEF-vdW^^ 

Yes 

0.06  +/-  0.03‘’  (6.0 

+/-  2.9) 

‘positive  value  of  the  difference  indicates  CH4  binds  more  strongly  on  the  Bronsted  site,  ’’energy 
values  are  non-ZPE  corrected,  ‘’estimated  Bayesian  statistical  error  calculated  from  a  distribution 
of  total  energies  estimated  using  an  ensemble  of  parameter  values  representing  the  functional. 
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Table  5;  Enthalpy  of  adsorption  of  CH4,  O2,  and  Ar  in  the  two  adsorption  regions  of  MOR  at  195 
K.  All  values  are  in  kJ/mol. 


Adsorbate 

12  MR  ehannel 

8  MR  side  poeket 

CH4 

-17.0 

-27.5 

O2 

-14.7 

-19.9 

Ar 

-14.0 

-20.6 

Table  6:  Enthalpy  of  adsorption  of  CH4,  O2,  and  Ar  in  ZSM-12  at  195  K.  All  energy  values  are  in 
kJ/mol. 


Adsorbate 

Binding  energy 

CH4 

-20.3 

O2 

-13.8 

Ar 

-16.6 

20 


(a) 


(b) 

Figure  1.  Differential  heats  for  CH4  adsorption  on  the  different  H-ZSM-5  samples  at  195  K.  (a) 
(•)H-ZSM-5(50)  and  (o)H-ZSM-5(280).  (b)  (♦)  ZSM-5(F,Si)  and  (0)H-ZSM-5(F,A1)- 
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Coverage  (mmol/g) 


Figure  2.  Differential  heats  of  adsorption  for  CH4  at  195  K  on  (■)  H-MOR  and  ( A)H-ZSM-12. 
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Figure  3.  Differential  heats  of  adsorption  for  O2  at  195  K  on  (■)  H-MOR  and  (A)  H-ZSM-12. 
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Figure  4.  Differential  heats  of  adsorption  for  Ar  at  195  K  on  (■)  Fl-MOR  and  (A)  H-ZSM-12. 
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Figure  5.  Most  stable  adsorption  structure  of:  (a)  methane  in  12  MR  channel,  (b)  methane  in 
8MR  side  pockets,  (c)  oxygen  in  12  MR  channel,  (d)  oxygen  in  8  MR  side  pockets,  (e)  argon  in 
12  MR  channel,  and  (f)  argon  in  8  MR  side  pockets  of  MOR.  Key:  Oxygen  (•),  silicon  (O), 
carbon  (•),  hydrogen  (O),  and  argon  (O). 
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Figure  6:  Most  stable  adsorption  structure  of  (a)  methane,  (b)  oxygen,  and  (c)  argon  in  the  12 
MR  channels  of  ZSM-12.  Key:  Oxygen  (•),  silicon  (O),  carbon  (#),  argon  (O),  and  hydrogen  (O 
). 


